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Magnetizations are introduced to the Generalized Random Energy Model (GREM) and nu- 
merical simulations on ac susceptibility is made for direct comparison with experiments in glassy 
materials. Prominent dynamical natures of spin glasses, i.e., memory effect and reinitialization, 
are reproduced well in the GREM. The existence of many layers causing continuous transitions 
is very important for the two natures. Results of experiments in other glassy materials such 
as polymers, supercooled glycerol and orientational glasses, which are contrast to those in spin 
glasses, are interpreted well by the Single-layer Random Energy Model. 
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§1. Introduction 

In spin glasses, it is well known that dynamical behavior strongly depends on history of the 
system after quenching from above the transition temperature T c . These phenomena are called 
aging and have been studied with various experimental protocols such as the isothermalJ3i) the 
T-shifti'i and the T-cycle one.BBQ'Era' From the theoretical point of view, aging phenomena have 
been studied along two different pictures, i.e., so-called droplet pictured 0' 00 and hierarchical 
pictureBIZm' The Generalized Random Energy Model (GREM)B@00) is a model that belongs 
to the latter. This model has a hierarchical structure causing continuous transitions as the system is 
cooled down. These continuous transitions correspond to successive branching process of free energy 
in the hierarchical picture. Bouchaud and Dean0 have shown that aging naturally occurs in this 
model and the time correlation function C(t + i w , t w ) satisfies a t/t w scaling law. Furthermore, we 
have recently made simulationdll0) similar to experiments on aging phenomena with temperature 
variations. As the consequence, it has been shown that results of experiments are reproduced well 
according to the hierarchical picture. 

Recently, Jonason et have made a new experiment, in which curious dynamical natures of 
spin glasses are observed very clearly. This experiment consists of the following two runs. In the 
first run, the sample is continuously cooled from T max (> T c ) to T m \ a (< T c ) at a constant rate, 



l 



and is immediately reheated at the same rate. During the cooling and the reheating, out-of-phase 
ac-susceptibility x" is measured as a function of temperature. We call this curve as x" e f- The 
difference between the observed value in the cooling and that in the reheating is not observed in 
spin glasses, while the hysteresis is observed in other glassy materials, such as polymer glasses, 
orientational glasses© and disordered ferromagnets@ In the second run, the sample is cooled 
from T max to a waiting temperature T wa i t (T m j n < T wa i t < T c ), and is kept at T wa ; t during a certain 
time interval. The sample ages and x" decreases during the interval. Then the system is cooled to 
T min and is reheated to T max without any stops. Hereafter this curve is denoted as Xwair 

An important result of this experiment is that Xwait merges with x" e f at very beginning of the 
resumed cooling as if the system forgets the aging at T wa jt (reinitialization of aging). But this 
aging is still imprinted and a dip of x'wa,it created in the cooling stage is exactly recovered in the 
reheating stage (memory effect). As Hammann et al have pointed out,© these two effects make 
dynamics in spin glasses temperature specific in the sense that the effect of aging at T wa ; t (or the 
difference between x" e { and Xwait) om y appears near T wa j t in the both cooling and reheating stages. 
The purpose of this manuscript is to study aging phenomena of the GREM with this protocol. 

The organization of this manuscript is as follows: In §2, the GREM is explained and magneti- 
zations are introduced to this model. In §3, the results of the simulations are presented. In §4, it 
is discussed how a variety of aging phenomena observed in glassy materials are interpreted within 
the GREM. 

§2. Model 

The GREM is schematically shown in Fig. 1. This model consists of L layers which are piled 
up hierarchically. The bottom points represent accessible states of the system and each branch 
represents a barrier over which the system goes to another state. Each branching point has N 
branches, so that the number of states is N L . It is assumed that ./V is large enough. Energy 
barriers of the n-th layer counted from the bottom, E n , are given randomly and independently 
according to the distribution 

(2.1) 



p n (E n ) = — ^-exp 



T c (n)_ 

where T c (n) is the transition temperature of the n-th layer. From the distribution, the averaged 
relaxation time of the n-th layer (r(n)) is easily calculated as 

(r(n))= dE n p(E n )T eME n /T) 
Jo 

oo (T< T c (n)), 

where to is a microscopic time scale. This means that a transition from the ergodic phase to the 
non-ergodic phase occurs at T c (n) in the n-th layer. The transition temperatures are chosen so as 
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to satisfy the inequality T c (l) < T c (2) < ••• < T C (L). Therefore, transitions occur continuously 
from the uppermost (the L-th) layer to the lowest one. 

Next, let us explain how magnetizations are introduced to the GREM . It is natural to assume 
that the nearer two states a and (3 are the stronger the correlation between the two magnetizations 
is, and that the distance between the two states d(a, (5) is measured by the layer from which they 
are separated, e.g., in Fig. 1 d(a,{3) = 1 and d(a,'y) = 2. To incorporate these aspects, we assign 
the value of the magnetization M a to state a as 

M a = M {a) + Mi(ai) + • ■ ■ + -Ml-iK-i), (2.3) 

where is the fe-th ancestor of a and A4k(ctk) is a contribution from the branching point. The 
value of A4k{oik) is given independently and randomly from distribution Dk{Mk) with the mean 



value Mk = 0. If d(a,/3) = k, the correlation between M a and Mg comes from the common 
contributions of M. n (n = k, k + 1, ■ • •) to these magnetizations 

L-l 

MjTp = Y,Ml. (2.4) 

n=k 

It decreases monotonically as k increases and the barrier between two states becomes higher, which 
is observed in the SK model© In the present simulation, the same uniform distribution is chosen 
for M. of all the layers as 

DniMn) J^ «*.ls*>. „ 5) 

1 (\M n \>^). 
The range of the distribution is chosen so that the variance of M a is independent of L. 

The dynamics of the system is described by the master equation for the probability P a (t) of 
finding the system at a state a at time t, 

±P a (t) = £ W a3 (t)P 8 (t) - J2 W Pa (t)P a (t), (2.6) 

where W a s(t) is the transition rate at t for going from (3 to a in a unit time. The uniform 
distribution is chosen for the initial condition, i.e., P a (0) = N~ L for each a. This means that the 
system is quenched from an infinitely high temperature. The transition rate W a n{t) is given as 



W a p{t) = r 1 Y, Jjk \ ex P 

k=d(a,/3) I 



cxp 



Ent\E n ((3) + H(t)Mf, 
T(t) 



T(t) 

(2.7) 

where H (t) is magnetic field and -Ex+i is hypothetical energy whose value is infinity. The factor in 
the braces on the right hand represents the probability that the system can be activated to the k-th 
layer but not to the k + l-th layer, and the factor 1/N k represents the probability that the system 
falls into a. The details of how the dynamics is simulated in the limit N — > oo are presented in 
ref.lil. 
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§3. Results of Simulations 

The simulation is done in the GREM with L = 2, T c (l) = 0.5 and T c (2) = 1.0. The amplitude and 
the period of the applying ac-field are 0.1 and IOOto, respectively. The system is cooled at the rate 
of 2.0 x 10 -5 T/to, and is immediately reheated at the same rate in the case of x'Li measurement. 
The cooling is intermitted at T wa it = 0.7 for 1.0 x 10 5 ro when Xwait i s measured. The waiting time 
is comparable to the sweeping time from T max to T min . The temperature T wait is set so as to satisfy 

T c (l) < T wait < T c (2). (3.1) 

In Fig. 2 there is obvious hysteresis, which does not disappear even in the case of slower rates 
(to 2.0 x 10~ 6 T/to) and longer periods (to 10 3 to). An important result of the simulation is reini- 
tialization, i.e., Xwait merges with x" ef in the early stage of the resumed cooling. 

In order to visualize the effect of the aging at 2~wait more clearly, the difference between Xref a nd 
Xwait ^ s shown in Fig. 3. If we regard the difference as a measure of the aging effect, we notice 
that the system behaves as if it quickly forgets the aging in the cooling stage (reinitialization) and 
again remembers when the system is heated back near T wa it (memory effect), which implies that 
dynamics of the GREM is temperature specific. 

To reveal the mechanism of these effects, x'o an d Xi evaluated from A4o and Ai\ (x" = Xo + Xi) 
are plotted in Fig. 4. The rapid increase of x" just after the resumed cooling is brought by x'o- As 
for the relaxation at T wa it, the first layer is quickly equilibrated and x'o almost retains a constant 
value because T wait > T c (l), while x'{ decreases since T wait < T c (2). In this sense, x'i is the affected 
part and x'o is the unaffected one. As the resumed cooling goes on, the affected part decreases 
and the unaffected part increases because the peak of x'o an d that of x'{ are near T c (l) and T c (2) 
respectively. As the result, Xwait merges with x" e i as ^ the system forgets the aging at T wa it- But 
Xi again recovers a large contribution to x" an d the system remembers the aging when the system 
is reheated to T wa j t . 

§4. Discussion 

As Hammann et al have pointed out, lit dynamics in glassy materials seems to be classified into 
two distinct types. In the first type, memory effect and reinitialization, which make the dynamics 
temperature specific, are observed. Spin glasses belong to this type as mentioned in §[|. In the other 
type, memory effect exists but reinitialization does not exist. Since the system is never reinitialized 
after quenching from above T c , the dynamics is cumulative. In this dynamics, the time evolution in 
the vicinity of T c is very important. Polymers supercooled glycerol© and the orientational 
glasses©© belong to this type. 

Now let us discuss how these two distinct types of dynamics observed in glassy materials are 
interpreted within the GREM. Time evolution of energy distribution P n (E n ,t), which is defined as 
the probability density that the system is found at time t in one of the states whose energy of the n-th 
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layer is E n , is very different according as T > T c (n) or T < T c (n)M> In both cases, P n (E n , t) tries to 
approach the equilibrium distribution proportional to exp{E n /T)p n {E n ) = exp({^ — T ^ }E n ) with 
increasing time. But the sign of the exponent is different in the two cases. Consequently, P n (E n ,t) 
quickly converges to the equilibrium distribution P^ q (E n ) = — T ^ } exp({y — T ^ }E n ) if 
T > T c (n), whereas P n (E n ,t) has a peak moving to higher energy level with increasing time if 
T < T c (n) (time evolution of P n (E n ,t) in the case of T < T c (n) is shown in Fig. 6, which is 
mentioned more closely later on). In this model, the peak position indicates the age of the layer. 
This means that the age of layers which satisfy T > T c (n) is kept to be since the peak of P^ q (E n ) 
is located at E n = 0. 

Although the aging process in each layer is cumulative, the dynamics of the GREM is temperature 
specific. Now let us discuss what happens in the case of L 3> 1. For a given temperature T < T C (L), 



there exists the n-th layer which satisfies T c (n — 1) < T < T c (n). As discussed in ref. 18, the crucial 
point is the fact that the layers below n are quickly equilibrated and do not contribute to slow 
dynamics, while those above n are almost quenched and they behave as if the time evolution stops. 
This means that the n-th layer is the activated one and mainly dominates slow dynamics of the 
system. This activated layer changes with temperature, that causes temperature specific dynamics 
of the GREM. Memory of aging at a given temperature is stored in the corresponding activated 
layer as a peak position of P n (E n ,t). The memory is preserved while the system is cooled since 
the layer becomes frozen, and is reinitialized if the system is heated up to a certain temperature 
T + AT > T c (n) because the peak of P n (E n ,t) created at T is destroyed. 

In the Single-layer Random Energy Model (GREM with L = 1), the dynamics is cumulative 
because changes of the activated layer do not occur. In Fig. 5, the relaxation of x" during positive 
T-cycle in this model is shown. The relaxation highly proceeds while the system is heated up, which 
is usually observed in cumulative systems. 0* On the other hand, relaxation of x" is reinitialized 
by positive T-cycle in the GREM. Ill* 

This understanding of aging phenomena observed in glassy materials is just the same as that 
obtained from studies of mean-field spin glasses]^ in which it has been concluded that dynamics 
of mean-field models with a one-step replica symmetry breaking solution is cumulative and that 
with a full replica symmetry breaking solution is temperature specific. 

In order to see how activated layers and frozen ones age, an example of time evolution of P n (E n ,t) 
is shown in Fig. 6. After quenching from an infinitely high temperature, the system is kept at 
T = 0.45. The number of layers is three and the transition temperatures are set at T c (l) = 0.6, 
T c (2) = 0.8 and T c (3) = 1.0, so that the 1-st layer is the activated one and the 2-nd and the 3-rd 
layers are frozen. We notice that the shifting speed of peak position of P n (E n ,t) in the frozen layers 
is slower than that in the activated layer. In fact, it can be proved analytically that in the case of 
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t ^> 1 and T c (n — 1) < T < T c (n), the energy distribution P n (E n ,t) satisfies the scaling, 

P k {E k ,t) = P k (E k - R k Tlogt) (k > n), (4.1) 
1 (k = n), 

Rk ^l TT T ^ (42) 

Rw) {k>n) - 

The validity is clear from the scaling plots shown in the insets of Fig. 6. The scaling implies 
that the age of the k-th. {k > n) frozen layer at t is nearly equal to that of the activated layer 
at t Rk , if age of each layer is measured by the peak position of the energy distribution. This 
frozen property becomes very important at experimental time scale t exp , which is considered to be 
10 13 < i cxp < 10 18 in ordinary spin glasses in unit of the microscopic time of the system. 
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FIGURE CAPTIONS 



Fig. 1 Structure of the GREM with L = 2 and N = 5. The open circles represent accessible states 
of the system. 

Fig. 2 Out-of-phase ac-susceptibility x" measured in the GREM with L = 2, T c (l) = 0.5 and 
T c (2) = 1.0. The amplitude and the period of the applying ac-field are 0.1 and IOOto, respectively. 
The system is cooled (and reheated) at the rate of 2.0 x 10~ 5 T/to. In the case of Xwait measurement, 
the cooling is intermitted at T wa i t = 0.7 for 1.0 x 10 5 ro. The solid line, open circles and full circles 
correspond to x" e f! Xwait (cooling) and x wa it (reheating), respectively. 

Fig. 3 Difference between Xr' e f and Xwait m the cooling and the reheating stages. The open circles 
and the full circles correspond to the cooling and the reheating data, respectively. 

Fig. 4 Out-of-phase ac-susceptibility x'o an d Xi evaluated from Mo and M\. The solid line, open 
circles and full circles correspond to x" c fi Xwait (cooling) and Xwait (reheating), respectively. 

Fig. 5 Effect of positive T-cycle in the Single-layer Random Energy model (the GREM with L = 1) 
with T c = 1.0. The ac-field with the peak amplitude 0.1 and the period IOOto is applied for the 
measurement of x" ■ After quenching from an infinitely high temperature, the system is kept at 
T = 0.5. Then a positive temperature perturbation AT = 0.35 is applied at t\ = 4 x 10 3 ro and is 
switched off at t\ + £2 = 1-4 x 10 4 ro. In the inset, £2 part of data is omitted and t\ and £3 parts are 
connected for comparison with the unperturbed data (dotted line). 

Fig. 6 Time evolution of P n (E n ,t) of the GREM with L = 3, T c (l) = 0.6, T c (2) = 0.8 and 
T c (3) = 1.0. After quenching from an infinitely high temperature, the system is kept at T = 0.45 
and P n (E n ,t) is measured at t = 10 4 - 5 , 10 5 - , 10 5 - 5 , . . . , 10 8 r (from left to right). In the inset, 
P n (E n ,t) is plotted as a function of E n — R n Tlogt, where i?i = 1.0, R2 = 0.75 and R3 = 0.42 (see 
text). 
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